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1.  Introduction 

The  diffuse  reflectance  peaks  occuring  at  *>685  nm 
which  were  reported  by  Morel  and  Prleur^  over  water  bod- 
1 es  bearl ng  phytoplankton  have  been  the  subject  of  specu¬ 
lation  over  the  past  few  years.  Enhanced  scattering  from 
chlorophyll  bearing  cells  has  been  observed  for  many 
years2»3,4.  Even  differential  scattering  of  circularly 
polarized  light  by  chloroplasts  has  been  observed^. 

There  seems  to  be  little  doubt  that  anomolous  dispersion 
provides  a  suitable  model  for  explaining  the  enhanced 
scattering  which  occurs  on  the  long  wavelength  side  of 
the  chlorophyll  absorption  peaks.  Gordon^  first  Invoked 
anofflolous  dispersion  to  explain  the  enhanced  reflectivity 
at  685  nm  but  later  retracted  this  hypothesis  in  favor  of 
chlorophyll  fluorescence^. 

The  one-dimensional  model  of  radiative  transfer  has 
been  effectively  used  by- several  researchers  1 to  give 
fairly  reliable  reflectivities  over  a  wide  range  of  ocean 
parameters.  It  Is  the  purpose  of  this  paper  to  extend 
that  model  to  account  for  chlorophyll  fluorescence. 


II.  Theory 


Let  us  consider  a  medium  which  contains  a  uniform 
distribution  of  particles  which  can  scatter,  absorb,  and 
emit  radiation.  We  will  also  assume  that  the  radiation 
Is  restricted  to  either  the  upward  (U)  or  downward  (0) 
d1  rect Ion  (see  Fig.  1).  Let  ( w/m^-nm)  denote  the 

spectral  Irradiance  at  wavelength  and  depth  z  In  the 
upward  (U)  or  downward  (D)  direction  In  a  medium  which  Is 
Infinitely  deep.  We  will  also  adopt  some  of  the  notation 
used  by  Gordon^  In  what  Is  to  follow.  It  Is  easy  to  show 
that  the  net  flux  absorbed  ,F/\(2)  (W/m^),  by  a  layer 
of  thickness  tz  for  wavelengths  In  the  spectral  band  of 
width  AX^  Is  q1 ven  by 


F.(z)  =  - -TZ -  AzAX 


(1) 


- ^ 

We  will  use  the  prime  superscript  to  denote  quantities 
which  are  functions  of  excitation  wavelength.  The  num¬ 
ber,  N/^(z),  of  photons  absorbed  per  unit  time  per  unit 
area  Is  then  given  by 

(2) 


/  F. 

Na(z)  .  " 


he 


where  h  Is  Planck's  constant  and  c  Is  the  speed  of 
light.  The  quantum  efficiency  n{X^Xp)  is  defined  to 
be  the  ratio  of  the  rate  at  which  photons  are  emitted 


from  Az  with  bandwidth  AAp  to  the  rate  at  which  they 
are  absorbed  by  chlorophyll  within  bandwidth  ax^  To  find 
the  fraction  of  photons  absorbed  by  chlorophyll  we  must 
multiply  eqn.  (2)  by  where  Is  the  total 

absorption  coefficient  and  is  the  chlorophyll 

Q 

absorption  coefficient.  Therefore  the  number  of  photons 
emitted  per  unit  time  per  unit  area  Is  given  by 


Np(z) 


n(x',Ap) 


<chl 

-V 
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(3) 


Now  let  Jf(z)  (w/m2-m-nm)  be  the  volume  fluorescent 
flux  at  depth  z  In  the  medium,  then 


Jp{z)AzAX  X 

-  he 


(4) 


Combining  eqns.  (1),  (2),  (3)  and  (4)  we  get 


Jp(z)AXp 


XAX 


dz 


(5) 


We  will  assume,  as  did  Gordon^,  that  the  fluorescent 
emission  band  can  be  represented  by  a  Gaussian  spectral 
distribution  centered  at  Xqp  ■  685  nm  with  a  •  10.64  nm 
and  therefore 


(6) 


Jp(z)  =  n(X.Ap)  jr  ^ 
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exp[  - 
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As  pointed  out  by  Gordon^,  Forster  and  Livingston^  have 
found  experimental 1y  that  the  quantum  efficiency  of  chlo¬ 
rophyll  a  (in  vitro)  is  independent  of  X^for  x'  <  Xqf 
and  falls  rapidly  beyond  this  value.  We  will  therefore 
assume  it  to  be  a  step  function  for  360  <  700nm. 

Using  eqn.  (6)  we  will  define 


'  ®’'P  t 

✓  2»a^ 


T  ^(^p) 

2,^ 
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Now  that  we  have  obtained  the  source  function  we  can  pro¬ 
ceed  to  the  derivation  of  the  differential  equations  gov¬ 
erning  the  upward  and  downward  fluorescent  irradiance. 

It  should  be  noted  that  the  irradiance  at  the  fluorescing 
wavelength  due  to  scattered  solar  flux  is  separable  from 
that  due  to  fluorescing  particles.  Kattawar  and  Plass^'^ 
already  solved  the  one  dimensional  equation  for  scattered 
solar  flux.  We  will  make  use  of  most  of  the  notation  In 
that  article  in  the  following  derivation. 


Referring  to  Fig.  1  let  us  consider  adding  an  in¬ 
finitesimal  layer  of  thickness  ^2  to  the  bottom  of  a 
homogeneous  medium  of  total  thickness  z.  We  will  also 
assume  tha-  the  fluorescence  source  function  Is  symmetric 
In  Its  emission  I.e.,  equal  amounts  upward  and  downward. 
Let  r(z)  and  t(2)  denote  the  reflection  and  transmission 
operators  respectively  for  a  layer  of  thickness  z.  Let 
HI^(z)  and  hD(z)  denote  the  fluorescent  irradlances 
emerging  from  the  top  and  bottom  of  the  layer  respec¬ 
tively.  Also  let  0j  denote  the  extinction  coefficient 
at  the  fluorescing  wavelength.  Considering  all  processes 
up  to  first  order  In  Az  we  get 

h‘^(z+Az)  »  H°(z){1-BjA2)  +  H°{z)  B  r(z)BjAz  + 
H*^(z)FBjAz  +  Jp(z)  Az  r(z) 

2 

where  F  +  B  =  Wq  the  single  scattering  albedo.  It 
should  be  noted  that  we  have  neglected  the  self  Induced 
fluorescent  term  i.e.,  at  685  nm  we  have  assumed  that  the 
excitation  at  this  wavelength  Is  caused  solely  by  the  ex¬ 
ternal  source  and  not  the  Internal  source.  This  assump¬ 
tion  is  clearly  Justified  due  to  the  fact  that  the  flu¬ 
orescing  irradiance  Is  small  compared  to  the  exciting 
irradiance.  This  equation  has  a  slnple  physical 


Interpretation  (see  Fig.  1).  The  first  term  on  the  right 
hand  side  of  eqn.  (8)  represents  the  radiation  emerging 
from  the  bottom  of  the  layer  of  thickness  z  which  passes 
through  the  sublayer  az  unscattered  The  second  term  re¬ 
presents  the  radiation  emerging  from  z  which  backscatters 
In  Az  and  then  is  reflected  from  the  overlying  layer  and 
then  emerges  from  az  unscattered.  The  third  term  Is 
radiation  emerging  from  z  which  forward  scatters  In  az 
and  emerges.  The  fourth  term  represents  the  fluorescent 
source  emitting  upwards  which  Is  reflected  from  the  over- 
lying  layer  and  then  emerges  through  Az  with  no  scatter¬ 
ings.  The  final  term  represents  the  downward  fluorescing 
radiation  from  the  layer  Az.  Now  passing  to  the  limit  as 
Az*0  we  get 

■IH°(r)  .  (ltr(2))  -  [(l-F)  -  p(2)B]0^h'’(2)(,) 

Using  the  same  notation  we  used  for  the  one  dimensional 
solution  In  reference  3  with  the  exception  that  X  has 
been  changed  to  v  In  this  paper  to  avoid  confusion  with 
wavelength,  we  have 

l-F«l-UQ-t>B»a-fB,  (10*) 

BCl-exp(-?eyZ)], 

'■(»)  •  yjTfjex-pr-fBx'O" 


(10b) 
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and 

p  -  ■  ’‘^lexPC-^ejZ/Z) 

It  should  be  noted  that  lim  hD(z)  =  0. 

z*- 

The  differential  equation  for  the  upward  Irradlance, 
H^(z),  exiting  the  top  of  the  layer  of  thickness  z  Is 
easily  found  to  be 


-  t(z)[Jp(r)/2  ♦  B6jH“(z)] 


(M) 


Using  eqns,  (10c),  (13).  and  (7)  In  eqn.  (14)  and  per¬ 
forming  the  integration  we  obtain  the  following  result 
subject  to  the  initial  condition  F^(0)  =  0  and  passing 
to  the  limit  as 
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(15) 


The  factor  appearing  in  eqn.  (15)  is  the 

downward  irradiance  just  above  the  ocean  surface.  It 
should  be  noted  that  in  order  to  interchange  the  order  of 
the  Integration  and  summation,  which  was  used  to  evaluate 
the  Integrals  over  the  variable  z,  the  power  series  ex¬ 
pansion  must  converge  uniformly. 


This  was  rather  easy  to  establish  for  the  Integrand  re¬ 
sulting  from  eqn.  (14).  We  can  reduce  eqn.  (15)  to  a 
rather  simple  form  as  follows.  First  for  typical  ocean 
models  fi/l'2“10"®  and  therefore  we  need  only  the 
n»0  term  in  the  summation.  Also  since  (B0Y/<Pi'f2)  <<1 
we  can  also  neglect  the  second  summation  term.  We  are 
therefore  left  with 


H^-) 


K(^p)C 

J — ? — 
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700 
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360 
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d  0 

♦prr 


dx 


(16) 


Also  over  the  entire  spectral  range  B'<<a^  then 


/  /  /  /  / 

C  "2a, 


Using  these  results  in  eqn.  (16)  and  recognizing  the  fact 
that  oPf  =  (l-«o)  Py  *  0  ^ we  get 

I  d 


K(Xp)  700  Hj(x')dX^ 

^^OF  360  (P^+P^) 

fl  a 


(17) 


which  Is  a  remarkably  simple  result 


III.  Calculation 


To  test  the  results  obtained  we  used  the  ocean  model 
given  in  Table  1.  The  chlorophyll  a  absorption  coeffi¬ 
cients  were  the  in  vivo  values  obtained  from  Morel  and 
Prieur^,  In  Table  II  we  present  the  results  of  our  cal¬ 
culations  for  the  quantum  efficiency  n  for  the  same  val¬ 
ues  of  chlorophyll  a  concentration  presented  by  Gordon^. 
We  should  note  that  the  reflectivity  Rs  at  685  nm  for 
zero  chlorophyll  concentration  computed  from  the  one  di¬ 
mensional  model  (see  ref.  10)  is  B/¥2  ^nd  for  the  model 
used  in  Table  I  give  0.095%.  The  extrapolated  value 
shown  by  Gordon^  using  the  data  of  Morel  and  Prieur^  was 
0.13%.  The  efficiency  values  in  parenthesis  in  Table  II 
used  Rs  «  0.13%  whereas  the  first  values  used  = 

0.095%.  In  Table  II  we  also  present  a  comparison  of  the 
emission  computed  by  using  eqn.  (15)  the  exact  result, 
compared  to  the  approximate  equation  (17).  As  can  be 
seen  the  approximate  result  is  accurate  to  less  than  1% 
for  the  model  presented.  The  range  of  the  Integral  over 
wavelength  was  taken  from  360  to  700  nm.  The  values  pre¬ 
sented  are  for  the  center  of  the  emission  line,  namely 
685  nm.  The  quantum  efficiencies  n  obtained  from  this 
one-dimensional  model  are  slightly  lower  than  those  ob¬ 
tained  by  Gordon  using  a  quasi-single  scattering  approxi¬ 
mation  with  full  three  dimensional  transfer;  however  some 
of  the  difference  could  be  due  to  the  use  of  different 


ocean  models. 

Using  eqn.  (18)  an  Interesting  upper  bound  can  be 
placed  on  the  quantity  R.-Rs*  ^  weighted 

mean  for  the  quantum  efficiency  n  =  0.47t  then  the  upper 
bound  win  be  "1.09%  and  this  result  Is  only  dependent  on 
n*  the  spectral  dependence  of  the  chlorophyll  a  absorp¬ 
tion  coefficient,  and  downward  irradiance. 

IV.  Conclusion 

We  have  obtained  an  exact  one-dimensional  solution 
to  the  equations  of  transfer  for  chlorophyll  fluorescence 
from  an  infinitely  deep  ocean  with  a  uniform  distribution 
of  fluorescing  bodies.  We  have  also  found  that  the  gen¬ 
eral  solution  could  be  reduced  to  a  very  simple  expres¬ 
sion  for  the  case  where  absorption  dominates  over  scat¬ 
tering.  This  allowed  us  to  establish  an  upper  bound  on 
the  difference  in  reflectivities  for  a  very  high  chlo¬ 
rophyll  concentration.  The  results  we  obtained  for  the 
quantum  efficiency  are  still  low  compared  to  the  average 
value  of  5%  obtained  by  Kleffer^^.  Gordon  also  found 
this  to  be  the  case  for  the  quasi-single  scattering  anal¬ 
ysis.  We  agree  with  Gordon's  conclusion  that  the  problem 
may  be  In  the  spectral  dependence  of  the  In  vivo  absorp¬ 
tion  coefficient  for  chlorophyll  a.  The  reason  being 
that  photon  cel  1 u1 ar  absorption  Is  not  all  converted  to 
photon  chlorophyll  absorption.  We  also  feel  that 


anomolous  dispersion  can  not  be  ruled  out  considering  the 
experimental  evidence  already  presented  for  It.  In  a 
future  paper  we  will  use  a  simple  dispersion  model  for 
the  refractive  index  and  starting  ^rom  this  calculate 
spectral  scattering  coefficients  In  the  region  of  the 
chlorophyll  absorption  bands.  This  will  introduce  a  con¬ 
centration  dependent  scattering  coefficient  which  may 
well  provide  enough  scattering  to  produce  the  desired 
peak  at  685  nra.  The  more  crucial  test  is  the  shape  of 
the  reflectivity  spectrum  about  the  peak.  The  Gaussian 
structure  gives  added  credence  to  the  fluorescense  hypo¬ 
thesis  since  anomolous  dispersion  scattering  does*  not 
usually  exhibit  this  structure. 
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Table  I.  Scattering  (65)1  absorption  (ffa),  and  total  extinction  co¬ 
efficients  (Bt)  the  water,  hydrosol,  and  chlorophyll  components  of 
the  ocean  model  as  a  function  of  wavelength.  is  the  downward  ir- 

radlance  Just  above  the  ocean  surface. 


Water 

(m-1) 

Hydrosol 

(m-1) 

Chi orophy 1 1  a 
(m" 1 / (mg-m3 ) ) 

Ba  3.80-2+ 

9.10-4 

1.00-2 

;  :o 

0.826 

Bs  1.20-2 

3.00-2 

- 

Bt  5.00-2 

3.09-2 

1.00-2 

Ba  2.20-2 

8.50-4 

1.20-2 

r  )0 

0.909 

Bs  9.40-3 

3.20-2 

- 

Bt  3.14-2 

3.28-2 

1.20-2 

)0 

1.20 

Ba  1.71-2 

7.90-4 

1.50-2 

Bs  7.60-3 

4.00-2 

- 

By  2.47-2 

4.08-2 

1.50-2 

Ba  1.53-2 

7.30-4 

2.00-2 

420 

1.51 

Bs  6.10-3 

8.60-2 

- 

By  2.14-2 

8.67-2 

2.00-2 

Ba  1.45-2 

6.90-4 

2.35-2 

*,40 

1.59 

Bs  4.90-3 

9.09-2 

- 

By  1.94-2 

9.16-2 

2.35-2 

Ba  1.56-2 

6.40-4 

2.50-2 

30 

1.84 

Bs  4.10-3 

8.70-2 

- 

By  1.97-2 

8.76-2 

2.50-2 

Ba  1.76-2 

6.10-4 

2.25-2 

30 

1.87 

Bs  3.40-3 

8.33-2 

- 

Bi  2.10-2 

8.39-2 

2.25-2 

♦  1  .ie 

notation  -N  means 

that  the  number 

Is  multiplied  by 

10-N 

r 


X  (nm) 

HDp{*</m2-nm) 

water 

hydrosol 

chlorophyll  a 

Sa 

2.57-2 

5.70-4 

1.80-2 

500 

1.76 

6S 

2.90-3 

8.00-2 

- 

0T 

2.86-2 

8.06-2 

1.80-2 

Ba 

4.77-2 

5.40-4 

1.60-2 

520 

1.66 

BS 

2.40-3 

7.67-2 

- 

bt 

5,01-2 

7.72-2 

1.60-2 

Ba 

5,58-2 

5.10-4 

1.25-2 

540 

1.62 

BS 

2.10-3 

7.40-2 

- 

Bt 

5.79-2 

7.45-2 

1.25-2 

Ba 

7.08-2 

4.90-4 

9.00-3 

560 

1.54 

bs 

1.80-3 

7.14-2 

- 

Bt 

7.26-2 

7.19-2 

9.00-3 

Ba 

1.08-1 

4.60-4 

7.50-3 

580 

1.55 

8$ 

1.60-3 

6.90-2 

- 

Bt 

1,10-1 

6.95-2 

7.50-3 

Ba 

2.44-1 

4.40-4 

7.00-3 

600 

1.51 

Bs 

1.40-3 

6.66-2 

- 

Bt 

2.45-1 

6.70-2 

7.50-3 

Ba 

3.09-1 

4.20-4 

7.50-3 

620 

1.47 

Bs 

1.20-3 

6.45-2 

- 

Bt 

2  45-1 

6.70-2 

7.50-3 

Ba 

3.29-1 

4.00-4 

1.00-2 

640 

1.44 

Bs 

1.00-3 

6.25-2 

m 

Bt 

3.30-1 

6.29-2 

1.00-2 

Ba 

4.00-1 

3.90-4 

1.50-2 

660 

1.40 

Bs 

8.00-4 

6.06-2 

m 

Bt 

4.01-1 

6.10-2 

1.50-2 

I 

I 

IIIIII . . 

. iH, 


Table  I  continued 


X  (nm) 

H0Q(w/m2.nm) 

water 

hydrosol 

chlorophyl 1 

680 

1.36 

Sa 

4.50-1 

3.70-4 

1.90-2 

3S 

7.00-4 

5.88-2 

- 

3T 

4.51-1 

5.92-2 

1.90-2 

700 

1.31 

3a 

6.50-1 

3.60-4 

2.00-3 

3S 

7.00-4 

5.72-2 

- 

3t 

6.51-1 

5.76-2 

2.00-3 

] 

1 

I 

I 

I 

1 


I 


Table  II.  Comparison  of  the  exact  emission  using  eqn.  (15)  with 
the  approximate  calculation  using  eqn.  (17). 


R  .  Rj 


chi  a(mg/m3) 


n  it) 


%  error 
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FIGURE  CAPTIONS 


Figure  1.  (a)  Schematic  representation  of  source  streams  giving 

rise  tc  fluorescent  radiation.  (b)  Schematic  repre¬ 
sentation  of  all  first  order  processes  encountered  upon 
addition  of  an  Infinitesimal  layer  of  thickness  Az  to  a 
layer  of  thickness  z  for  the  downward  Irradlance. 


